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On the basis of a general formula obtained earlier via fourth-order perturbation theory within 
the framework of macroscopic quantum electrodynamics, the van der Waals potential between two 
neutral, unpolarized, ground-state atoms in the presence of a homogeneous, dispersing and absorbing 
magneto-electric sphere is studied. When the radius of the sphere becomes sufficiently large, the 
result asymptotically agrees with that for two atoms near a planar interface. In the opposite limit 
of a very small sphere, the sphere can effectively be regarded as being a third ground-state atom, 
and the nonadditive three-atom van der Waals potential is recovered. To illustrate the effect of 
a sphere of arbitrary radius, numerical results are presented for the triangular arrangement where 
the atoms are at equidistance from the sphere, and for the linear arrangement where the atoms 
and the sphere are aligned along a straight line. As demonstrated, the enhancement or reduction 
of the interaction potential in the presence of purely electric or magnetic spheres can be physically 
understood in terms of image charges. 

PACS numbers: 12.20.-m, 42.50.-p, 34.20.-b, 42.50.Nn 



I. INTRODUCTION 



Electromagnetic and material ground-state fluctua- 
tions are well known to produce observable effects such 
as dispersion forces between atoms, between atoms and 
bodies, or between bodies [1]. The van der Waals (vdW) 
potential between two ground-state atoms in free space 
was first calculated by London for small interatomic sep- 
aration (nonretarded limit), using perturbation theory to 
the leading (second) order [2]. In this limit, the result is 
an attractive potential proportional to r~ 6 , where r de- 
notes the interatomic distance. The (nonretarded) Lon- 
don potential was extended to arbitrary distances be- 
tween the two atoms by Casimir and Polder [3], using 
fourth-order perturbation theory within the framework 
of full quantum electrodynamics (QED). In particular for 
large interatomic separation the potential was predicted 
to vary as r~ 7 . Recently, the closely related Casimir in- 
teraction between two magnetoelectric spheres has been 
studied by means of a scattering method [4], where the 
inclusion of higher-order multipoles have been shown to 
lead to corrections of the Casimir-Polder result. In the 
three-atom case, a non-additive term prevents the po- 
tential from just being the sum of three pairwise contri- 
butions. This three-atom vdW potential was calculated 
first in the nonretarded limit by pursuing the perturba- 
tion calculation to the third order [5], and then for ar- 
bitrary interatomic distances, by using sixth-order per- 
turbation calculation [6]. Later on, a general formula for 
the non-additive A-atom vdW potential was derived by 



summing up the response of each atom to the quantized 
field caused by the other atoms [7] and by calculating the 
difference in the zero-point energy of the electromagnetic 
field of a large cavity with and without the atoms [8] . 

The presence of macroscopic bodies modifies the fluc- 
tuation of the electromagnetic field, and consequently, 
the interatomic vdW interaction. A general formula ex- 
pressing the vdW potential between two ground-state 
atoms in the presence of electrically polarizable bodies 
in terms of the Green tensor of the body-assisted elec- 
tromagnetic field was first obtained by means of lin- 
ear response theory [9], an later by treating the effect 
of the bodies semiclassically [10]. An exact derivation 
of the formula based on fourth-order perturbation the- 
ory within the framework of macroscopic QED has been 
given recently [11], where both electrically and magneti- 
cally polarizable bodies — referred to as magneto-electric 
bodies — are explicitly taken into account, and a general- 
ization of the formula to the A^-atom case has also been 
given [12]. The two-atom potential in the presence of ar- 
bitrary magnetoelectric bodies has recently been gener- 
alized to atoms having both electric and magnetic polar- 
izabilities using linear response theory and was explicitly 
calculated for atoms embedded in a bulk medium [13]. 
Various special cases such as two atoms placed between 
two perfectly conducting plates [10], in a bulk magneto- 
electric medium [11, 14-16], near a perfectly reflecting 
plate [11], and near a planar magneto-electric multilayer 
system [11] have been considered. Local-field corrections 
to the vdW potential that appear when the atoms are 
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embedded in optically dense media have also been ad- 
dressed [17, 18]. 

On the experimental side, the vdW interaction be- 
tween a single atom and a body has been explored by 
means of detecting the intensity of an atomic beam trans- 
mitted through a parallel-plate cavity [19]; direct force 
measurement using atomic-mirror techniques [20, 21]: 
measuring the intensity of a diffracted atomic beam from 
a transmission-grating [22] ; making use of quantum re- 
flection from a solid surface at nonretarded [23] and re- 
tarded [24, 25] atom-surface separations; determining 
their effect on the collective oscillation frequency of the 
magnetically trapped atoms [26]. Observations of inter- 
atomic vdW interactions based on a determination of the 
scattering cross sections in the atomic collisions between 
two ground-state atomic beams [27], between an atomic 
beam and the atoms of a stationary target gas [28, 29], 
and between a beam of ground-state atoms and a beam 
of excited atoms [30] have been reported. 

As already mentioned, the theoretical studies of 
medium-assisted interatomic vdW interactions have so 
far concentrated on bulk media and infinitely extended 
planar (multi-layer) bodies Here we shall consider the 
vdW interaction between two ground-state atoms located 
near a finite-size body, namely, a sphere. With recent 
progress in fabrication of metamaterials in mind, we al- 
low for the sphere to exhibit both electric and magnetic 
properties. Note that the single-atom vdW potential in 
the presence of a sphere has been investigated earlier [31]. 

The paper is organized as follows. The basic formulas 
for calculating the vdW potential between two atoms in 
the presence of an arbitrary arrangement of magneto- 
electric bodies are summarized in Sec. II. In Sec. Ill, 
the theory is applied to the case of a magneto-electric 
sphere, and the limiting cases of large and small sphere 
are considered. Detailed numerical results are presented 
in Sec. IV. Finally, the paper ends with a summary in 
Sec. V. 



II. BASIC FORMULAS 

Consider two neutral, unpolarized, ground state atoms 
A and B with spherically symmetric polarizabilities in 
the presence of an arbitrary arrangement of dispersing 
and absorbing magneto-electric bodies. The total force 
acting on the atoms can be derived from the potential 

U{v a ,t b ) = Ua{va) + U b {vb) + Uab{va,vb) (1) 
according to 

F^(r A ,) = -V rA ,£/(r A ,r s ), (2) 
(A' = A, B), where U A > is the single-atom potential [31] 

U A >(r A ,) = / duu 2 a A >(iu)TrG {1) {v AI ,Y A ,,iu), 
2tt Jo 



and Uab is the two-atom interaction potential [11] 
hfi 1 f°° 

U AB (r A ,r B ) = — — / duu 4 a A (iu)a B (iu) 
27T Jo 

x Tr[G(T A ,r B ,iu)-G(TB,TA,iu)]. (4) 

In Eqs. (3) and (4), a a> is the (lowest-order) polarizabil- 
ity of the atom A' 
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a A > (uj) = lim — V . - 

rj^0+ 3ft ^ (UJ A ,) 2 - UJ 2 - 17]U! 



(5) 



with u A °, and d A °, being, respectively, the transition fre- 
quency and transition electric dipole moment between 
the kth excited state and the ground state of atom A 1 , 
and G(r, r',u) is the classical Green tensor obeying the 
differential equation 



V x 



1 



V x G{r,r',u) 



-— £ (r,u;)G(r,r>) = 5(r-r') (6) 

together with the boundary condition at infinity. Note 
that the electromagnetic and geometric properties of the 
bodies are fully incorporated in the Green tensor via the 
space- and frequency-dependent permittivity e(r, u)) and 
permeability /i(r, to). When r and r' denote two positions 
in free space which can be connected without crossing a 
body, then the Green tensor can be decomposed as 

G(r, r', W ) = G<°> (r, r', w ) + G« (r, r', w ), (7) 

where G^ is the free-space Green tensor which is ob- 
tained from Eq. (6) by letting e(r,w) = ^{r,uj) = 1, and 
G^ is the scattering part of the Green tensor. 

In what follows we concentrate on the atom-atom 
interaction potential U AB in the presence of a sphere 
(for the single- atom potentials U A (r A ) and U B (r B ), 
see Ref. [31]). According to Eq. (7), the potential 
Uab(ta,*b) can be cast in the form 

U AB {v A ,v B ) = U^{v A ,v B ) + U {h \v A ,v B ), (8) 

where 

roo 

UW(r A ,r B ) = -^p / duu^aAii^asiiu) 
An Ja 

x Tr[G^(r A ,r B ,iu) ■ G^(r B , r A , iu)] (9) 

is the potential observed in the case when the two atoms 
are in free space, and is the body-induced part which 
can be written as 



(3) 



U^\r A ,r B ) = f7 (1) (r A ,r B ) + {^A,r B ), 
where 

^2 roc 

U^\r A ,r B ) = / duu 4 a A (iu)a B (iu) 

n Jo 

x Tr[G (0) (r A ,r B ,iu) ■ G { ~ 1 \y b y Ai iu) 



(10) 



(11) 
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is the contribution due to the cross term of the free-space 
part and the scattering parts of the Green tensor, and 

ha 2 f°° 

U (2) (r A ,r B ) = — / duu 4 a A (iu)a B (iu) 
27T Jo 

x Tr[G^(v A ,r B ,iu) ■ G^{r B ,r A ,iu)\ (12) 

is the scattering-part contribution. 

The free-space Green tensor reads (see, e.g., Rcf. [32]) 



G (a \r A ,r B ,Lu 



4iru 2 P 



f(-ilw/c)I - 



,ilu>/c 



(13) 



(r A ^r B ), where I denotes the unit tensor, 1 = r B —r A , 
I = 111, and 



f(x) = l + x + x 2 , 
g(x) — 3 + 3x + x 2 . 



(14) 
(15) 



By substituting Eq. (13) together with Eqs. (14) and (15) 
into Eq. (9), the well known Casimir-Polder interaction 
potential between two ground-state atoms in free space 
are obtained, 
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FIG. 1: Two atoms A and B in the presence of a sphere 

Hankel functions of the first kind, hn\x), and Legendre 
functions, P™(x), as follows: 

M nro ,±i(r,g) = T-7—7; h£\qr)P™ (cos 9) Sm (m^)eg 

blU u COS 

,(,), . dPI™(cos#) cos, ,, 
- h& (qr) — ^ 1 (m<P)et, (18) 



U^(r A ,v B ) = - 



16n 3 e 2 l e Ju 
6lu 5l 2 u 2 
c c z 



dua A (iu)a B (iu)e 2lu l c 
■ (16) 



2/ 3 M 3 l 4 u 4 



Needless to say that the scattering part of the Green 
tensor depends on the specific arrangement of the bodies 
under consideration. 



III. TWO ATOMS IN THE PRESENCE OF A 
MAGNETO-ELECTRIC SPHERE 

Let us consider two atoms A and B in the presence of 
a homogeneous magneto-electric sphere of permittivity 
e(lo), permeability and radius R. Choosing the 

coordinate system such that its origin coincides with the 
center of the sphere, we may represent the scattering part 
of the Green tensor as [33] 



G (1) {v A ,r B ,Lu) = — V \^ (n-mj. ^ ^ where j n (z) is the spherical Besscl function of the first 

4ttc ^-f n(n+l) ^ n (n+m)\ m v^a ^ = hH (h. = ,.,IA = . /w^W^W = «.rfV,iW anH 



tvt i \ n ( n + X ) 7,(1 
N„ m ,±i(r, q) = — '- h)l 

qr 



\qr)P^ (cos ef2(m<P)e r 



1 d[rh^(qr)} 



+ 



qr 



dr 



dP™(cos<9) cos 



sm 



(m<j))eo 



T-^P„ m (cos^W>)e 



sin( 



cos 



(19) 



with e r , eg, and being the mutually orthogonal unit 
vectors pointing in the directions of radial distance r, 
polar angle 9, and azimuthal angle <j), respectively (inset 
in Fig. 1). The coefficients B™(w) and B%(u) in Eq. (17) 
read 



H(u!)[zoj n (zo)}'jn(zi)-[zij n (zi)]'j„(zo) 



^l(uj)[z h^ ) (zo)}'jn(zi)~[zij n (z 1 )]'h ( n ) (z a ) 



(20) 



£(u)[z jn(zo)]'jn(zi) ~ [z 1 j n (z 1 )}' j n (z a ) 



e(aj)[z h ( n ) (zo)}'j n (z 1 ) - [z 1 j n (z 1 )}'h ( n\z ) 



(21) 



n(n + l) ^-^ (n+m) 



(r A ,uj/c)M nnitP (r B ,uj/c) 

p=±i 

+B»(u,)N (r A ,uj/c)N nm)P (r B ,(j/c)] , (17) 

where M nTO>p (r, q) and N nTO;P (r, q) are even (p=+l) and 
odd (p= — 1) spherical wave vector functions. In spherical 
coordinates, they can be expressed in terms of spherical 



kind, zo = kR (k = u>/c), z\ = e(uj)ii(lo)zq = n(tu)zo, and 
the prime denotes differentiation with respect to the re- 
spective argument. Without loss of generality, we assume 
that the two atoms are located in the xz-plane (Fig. 1), 



r A = (r A ,6 A ,0), r B = (r B ,6 B ,n). 



(22) 



Substituting Eqs. (18) and (19) [together with Eq. (22)] 
into Eq. (17), and performing the summations over m and 
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p, we derive (Appendix A) 

G^(v A ,r B ,iu) =J2G\f(r A ,r B ,iu)e iA e jB (23) 

ij 

(i, j = r, 8, </>), with the nonzero elements being 



Gi])(r A ,r B ,uj) 



ic 



^n(n+l)(2n + l)^HP n ( 7 )gl 1 ), (24) 



n=l 



2(2n+l)B^(w)i^(7)_( 3 ), (25) 



A.-kujta'tb 



n=l 



2n+ 1 



Attc ^-^ n(n + 1) 

n—l v ' 

p M ,/:(7) fi (i) | c 2 B»(u) F n ( 7 ) (4) 

B n i") Fn{l) QV + ^p^/i' 



where = #a + 0g, 7 = cos 0, and 

QW=hW(kr A )h£\kr B ), 



^=h^(kr A )[zh^(z)}' z= 



kr B ' 



(26) 



(27) 



(28) 



Qi 3 ^^^)^ 1 )^)];^, (29) 



F„( S )=n(n+l)P n ( S )-<(4 (31) 

To facilitate further evaluations, it is convenient to repre- 
sent the free-space Green tensor (13) in the same spher- 
ical coordinate system as the scattering part, so that its 
nonzero elements read 



G^}(r A ,r B ,iu) 



4iru 2 l 



G ( °J {8r) (r A ,r B ,iu) 



B [Z 2 /(0 cos6- 5 (0 l A l B ]e- 6 , 

(32) 

1 sin 6 



4ttu 2 1 5 

x [l 2 f^)±g{0r A{B) l A(B) ]e-^ (33) 



GfJ (r A ,r B , iu) 



4iru 2 l 5 

x [Z 2 /(0 cos9 - g(Z)r A r B sin 2 0] e"«, (34) 



G W (rA)rB>itt) = __^_ /(0c -C | (35) 

where £ = lu/c, l A — r B cos 6 — and = — cos e. 

Recalling Eqs. (10)— (12), we may write the body- 
induced part of the interaction potential as 



U {h \v A ,v B ) 
where 



^[u^(rA,r B ) + U^(v A ,r B )], (36) 



hy 2 r°° 



Gf(r A ,r B , iu)G\;> (r A , r B ,iu), (37) 



du u A a A {iu)a B {iu) 



.(i) 



fc 1 f 00 

U^\r Al r B ) = — ^-jT J duu 4 a A {iu)a B (iu) 



GV{T A ,T B ,iu) ' , (38) 



with G\^ (r A ,r B ,iu) and G^ (r^, r Bl iu) according to 
Eqs. (24)-(26) and Eqs. (32)-(35), respectively (see Ap- 
pendix B). Further evaluation of (ta^b) requires 
numerical methods in general. Before doing so, let us 
consider the limiting cases of large and small spheres. 



A. Large sphere 

The limiting case of a large sphere may be defined by 
the requirement that 



5 a> =r A > - R<.R (A' = A, B) 



(39) 



= \vh^{y% =k r A \^(z)]' z=krB , (30) and 



z<i? -» e = e A + e B < 1 (40) 

(cf. Fig. 1). In this limit, Eq. (36) leads to (Appendix C) 

U^(v A ,r B ) = 

I ^{[4X 4 -2^ 2 + X 2 ( , 2 + , 2 ) ]7 01 

- 3 £ (/° 2 + 4^ 2 )}< <«> 

where A = -i?e, (5 ± =6 B ±6 A , Z+ = (A 2 + 8\ ) 1 / 2 , and 

/>0O 

= / du a A (iu)a B (iu) 
Jo 

e(iu) — l" 1 



hi 



[26 + e(iu)fi(iu)R 2 u 2 /c 2 



e(iu) + 1 



• (42) 
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Note that the case of a purely electrically polarizable 
sphere can be simply obtained by setting /j,(iu) = 1 in 
Eq. (41). 

For a purely magnetically polarizable sphere [e(iu) = 1] 
Eq. (36) leads, under the conditions (39) and (40), to 
(Appendix C) 



U^(v A ,v B ) = 



h[6 2 _ -2X 2 + 3S+(l+ -5+)] 



Q4n 3 e 2 c 2 l 5 l + 

f°° , 2 , . . , \p(iu) - l\Miu) - 3] 

x / &uu 2 a A (iu)a B {iu y K ' - \ LPV - ' -. (43) 

Jo 



fi(iu) + 1 



As expected, Eqs. (41) and (43) are in agreement with 
those found for a half-space (Z + /i?— >0) [11]. 



B. Small spheres 

In the opposite limit of a small sphere, where 

i?«r A , (A'=A,B), (44) 
Eq. (36) leads to (Appendix C) 



U (b) (r A ,r B ) 

h 



647r 4 £fjr|r^ 3 Jo 



f-OO 

/ dua A (iu)a B (iu)e- (rA+rB+l '> u/c 
Jo 



x \ a sp (iu) 



m{g(b)[2(l+a)-g(a) sin 2 6] +2a 2 f(b)} 



+ 



{[(2l 2 ~r A r B cosG)f(a),f(b) + 2a 2 f(b)r A l / 



-2b 2 f(a)r B l B ] sin 2 <d-4(l + a)(l + b)l A l B cos 6} 



+ ^(l+o)(l + 6)A, p (i« 



(a = r A u/c, b = r B u/c), where 



f ^r A r B . 2 
-2/(0 cos6 



/ \ q e(^) — 1 
a sp {u) = 4ire R 3 — 



e(w) + 2' 

47ri? 3 - 1 
yu /u(o>) + 2' 



(45) 



(46) 



(47) 



Let us consider a sphere to which the Clausius-Mossotti 
relation applies, so that 



e(uj) + 2 3e 



(48) 



with n k and a k {ijj), respectively, being the number den- 
sity and the polarizability of the atoms of type k. In this 
case, Eq. (46) can be rewritten as 



= ^n k a k {uj) 



(49) 



where N k is the number of atoms of type k of the sphere. 
Accordingly, the magnetic analog of Eq. (49) is 



(50) 



with (3 k (cu) being the magnetizability of the atoms of 
type k. Hence, we may replace in Eq. (45) the sphere 
parameters a sp (iu) and f3 sp (iu), respectively, with the 
electric and magnetic polarizability of a single atom [say 
a k (iu) and j3 k (iu)}, to obtain the nonadditive interac- 
tion potential of three atoms, two of which being purely 
electrically polarizable whereas the third atom being si- 
multaneously electrically and magnetically polarizable. 
Indeed, after a straightforward but lengthy calculation, 
it can be shown that in the case of a purely electrically 
polarizable sphere, Eq. (45) [[3 sp (iu) = 0] leads to the in- 
teraction potential between three electrically polarizable 
atoms, as derived in Refs. [5, 7, 8, 34]. 

In the retarded limit where l,r A ,r B 3> c/cj m ; n (w m ; n 
denoting the minimum frequency among the relevant 
atomic and medium transition frequencies), due to the 
presence of the exponential term in the integral in 
Eq. (45) , only small values of u significantly contribute to 
the integral. Therefore, the electric and magnetic polar- 
izabilities can be approximately replaced with their static 
values. After performing the remaining integral and ex- 
pressing all the geometric parameters in terms of r A , r B , 
and I, we arrive at 



U {b \r Al v B ) 



hca A (0)a B (0) 



32w i £ 3 r A r B l 5 (r A + r B + l) 7 
a sp (0){5[ft,i(rA,r B ,Z)] + S[h\(r B , I, r A )\ 

2 2 

+ S[h 1 (l,r A ,r B )]} + 1 ^-(3 sp (0)S[h 2 (r A ,r B ,l)} 



(51) 



where 



hi(x,y,z) = 3x 6 y 2 {y-x){x + y+7z){x 2 + 7xy+lly 2 ) 
- x 4 y 2 z 2 {53x 4 + 280x 3 y - Yi7x 2 y 2 

- 329xy 3 - 623xy 2 z - 192y 2 z 2 ), (52) 

h 2 (x,y,z) = 3x 4 (y-x)(x + y+7z)(x 2 + 7xy + lly 2 ) 
-2x 3 z 2 (x + y){26x 2 + 93xy-l33y 2 )-7x 2 z 5 (3x-2y) 
- 14x 3 z 3 (2x 2 - 3xy-Uy 2 ) -x 3 z 4 (l7x+miy) 

+ 2xz 6 (31a; + 105y) + 5z 7 (14a; + z), (53) 

and S[f(x, y, z)] = f(x, y, z) + f(y, x, z). 

In the nonretarded limit where l,r A ,r B <Cc/[n(0)aj max ] 
(k-Wx denoting the maximum frequency among the rel- 
evant atomic and medium transition frequencies), the 
leading contribution to the integral in Eq. (45) comes 



6 



from the region where e (r A +r B +l)u/c ^ ^ so Eq (45) 
reduces to 



U {b Hr A ,r B ) 



647r 4 egr^r|P 



1 



1 - 72 (4W b +r^rs siri 9) cos 9 +cos^ 9 



+ 



Ji 



r A r B ( r A r B • 2 2 o i / 
— j- I -7^ sm 9--cos9 ) J 2 



where 



Jo 

-/ 

Jo 



du a a (iu)as (iu)a sp (iu) , 



du u 2 aA{iu)aB{iu)[} sv {iu). 



(54) 



(55) 



(56) 



In particular, in the case of a purely electrically polariz- 
able sphere (J 2 = 0), Eq. (54) can be written in a very 
symmetric form. For this purpose we introduce the unit 
vectors a, b, and c pointing in the directions of ya, 1, and 
— tb, respectively (see Fig. 2). Noting that l A and Ib de- 
fined below Eq. (35) are the components of the vector 1 
in the directions of ta and yb and can thus be written 
as l(a-b) and —l(b-c), respectively, we see that 

- ^ (<UaIb + r A r B sin 2 9) +cos 9 = 3 (a ■ b) (b • c) (57) 
and can rewrite Eq. (54) as 

^^' r ^= 64^P [l-3(a.6 ) (6.c ) (ca ) ]J 1 . 

(58) 

If a sp (iu) in Eq. (55) is again identified with the the elec- 
tric polarizability of a single atom, Eq. (58) is nothing but 
the formula for the nonretarded three-atom interaction 
potential, which was first given by Axilrod and Teller [5] . 



IV. NUMERICAL RESULTS 

The effect of a medium-sized magneto-electric sphere 
on the mutual vdW interaction of two identical two-level 
atoms is illustrated in Figs. 3 and 4 showing the ratio 
Uab /t/ (0) [cf. Eq. (8)]. The results have been found 
by exact numerical evaluation of Eq. (36) together with 
Eqs. (B1)-(B7), where the permittivity and permeability 
of the sphere have been described by single-resonance 
Drude-Lorentz models, 

s(uj) = 1 + - , (59) 

^) = 1 + - . (60) 

In Fig. 3, a configuration is considered where the two 
atoms are positioned at equal distances from the sphere, 




FIG. 2: The triangle formed by the two atoms (at corners A 
and B) and the sphere (at the corner C) in the small-sphere 
limit. It is seen that the vector products a b, b e and c a 
in the Axilrod and Teller's formula [5] are equal to — cos a, 
— cos/3, and — cosO, respectively. 



ta = tb, briefly referred to as triangular configuration, 
and Uab /U^ is shown as a function of the angular sep- 
aration 9 for three different values of the atom-sphere 
separation. For a purely electrically polarizable sphere 
[Fig. 3(a)], depending on the separation angle between 
the atoms, a (compared to the free-space case) relative 
reduction or enhancement of the vdW potential is pos- 
sible, while for a purely magnetically polarizable sphere 
[Fig. 3(b)], the potential is typically reduced [note that 
for very small angular separations, a slight enhancement 
is possible, as can be seen from the inset in Fig. 3(b)], 
and the reduction increases with the angular separa- 
tion. In both cases, the sphere-induced modification is 
strongest when the atoms are at opposite sides of the 
sphere (9 = n). Note that for small atom-sphere sepa- 
rations (solid curves) and small angular separations, the 
potential qualitatively agrees with the potential obtained 
for two atoms placed in parallel alignment near a semi- 
infinite half space [11], as expected from the results in 
Sec. Ill A. 

In Fig. 4, a configuration is considered where the two 
atoms and the sphere center are aligned on a straight line, 
briefly referred to as linear configuration, and Uab/U^ 
is shown as a function of the interatomic distance for 
three different values of the position ta of atom A which 
is positioned between the sphere and atom B. Unless 
both atoms are very close to the sphere, the sphere gives 
always rise to a (compared to the free-space case) relative 
enhancement of the vdW potential between the atoms; 
only for very small atom-sphere separations the poten- 
tial can be reduced if the sphere is purely magnetically 
polarizable [inset in Fig. 4(b)]. Figure 4(a) shows that in 
the presence of a purely electrically polarizable sphere the 
relative enhancement of the potential increases with the 
interatomic separation I and approaches a limit for larger 
interatomic separations, which depends on the separation 
distance between atom A and the sphere. From Fig. 4(b) 
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FIG. 3: The mutual vdW potential of two identical two- 
level atoms in a triangular configuration with (a) a purely 
electrically polarizable sphere with a>p e /o>io =3, WTe/^io = 1, 
and 7e/u>io = 0.001 and (b) a purely magnetically polarizable 
sphere with wp m /o>io =3, WT m /wio = 1, and 7 m /o»io = 0.001 
is shown as a function of the atom-atom angular separation 
G (wio is the atomic transition frequency). The sphere ra- 
dius is R = c/ojio and the distances between the atoms and 
the center of the sphere are va = vb = 1.03c/wio (solid line), 
1.3c/wio (dashed line), and 2c/cjio (dotted line). is the 
potential observed in free space. 




(a) 



0.01 0.1 1 10 100 1000 




0.01 0.1 1 10 100 1000 

lijJw/c 



FIG. 4: The mutual vdW potential of two identical two-level 
atoms in a linear configuration with (a) a purely electrically 
polarizable sphere and (b) a purely magnetically polarizable 
sphere is shown as a function of the interatomic distance I. 
Atom A is held at a fixed position between atom B and 
the sphere center with ta = 1.03c/wio (solid line), l.lc/wio 
(dashed line), and 1.3c/wio (dotted line). All other parame- 
ters are the same as in Fig. 3. 



it is seen that in the presence of a purely magnetically 
polarizable sphere the relative enhancement of the poten- 
tial increases with the interatomic separation I, reaches 
a maximum, and decreases with a further increase of I. 
In agreement with the results of Sec. Ill A, the poten- 
tial observed for small atom sphere separations (solid 
curves) and small interatomic separations qualitatively 
agrees with the potential obtained for two atoms placed 
in vertical alignment near a semi- infinite half space [11]. 

Many features of the vdW potential observed in Figs. 3 
and 4 can be subject to a physical interpretation via the 
method of image charges (the same approach as has been 
used for a planar geometry Ref. [11]). Although being 
strictly valid only for sufficiently small atom-atom and 
atom-surface distances (such that retardation is negligi- 
ble) and being most easily applicable in the perfect con- 
ductor limit, this approach yields qualititive predictions 



for the sphere- induced enhancement and reduction of the 
potential which apply beyond this case. According to the 
image-charge method, the effect of the boundaries is sim- 
ulated by suitably placed image charges of appropriate 
magnitudes, so that the two- atom vdW potential effec- 
tively consists of interactions between fluctuating dipolcs 
A and B and their images A' and B' in the sphere, with 

Hint - V AB + Vab> + V BA , (61) 



being the corresponding interaction Hamiltonian. Here, 
Vab denotes the direct interaction between dipole A and 
dipole B, while Vab 1 and Vba 1 denote the indirect inter- 
action between each dipole and the image induced by the 
other one in the sphere. The leading contribution to the 
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FIG. 5: Two electric dipoles near a purely electrically polar- 
izable sphere are shown in a triangulare configuration where 
the dipole-dipole angular separation (a) is small enough to 
replace the sphere, approximately, by a halfspace and (b) is 
equal to 7T. 




FIG. 6: Two magnetic dipoles near a purely electrically po- 
larizable sphere are shown in a triangulare configuration for 
a (a) very small angular separation (b) maximal angular sep- 
aration. 



energy shift is of second order in H- ln t, 
AE AB - -2, fiK+ ^) 

n,m D 

X (n A ,TO B |#int|0A,0i3), (62) 

where In^m)) denotes the energy eigenstates of atom 
A(B) with eigenenergies <^a(b) and the prime indicates 
that the terms Ua(b) = ®A(B) arc n °t included in the 
sum. 

In this approach, corresponds to the product of 
two direct interactions and is negative in accordance with 
Eq. (62). is the product of two indirect interactions 
and is also negative. The terms containing one direct 
and one indirect interaction are contained in Since 
the total potential is equal to U A b = + U {1) + U {2 \ 
the sum + represents the effects of the medium. 
The relative signs and strengths of f/W and will 
determine whether the free space vdW interaction is en- 
hanced or suppressed. The orientations of the dipoles 
A and B are random and independent of each other, so 
that strictly speaking the signs of all dipole-dipole in- 
teractions has to be obtained by averaging over all pos- 
sible orientations. The effect of such averaging on the 
sign of the interactions can be reproduced by restricting 
the attention to the maximally attractive case of both 
dipoles pointing in the same direction parallel to their 
connecting line, with the dipole-dipole interaction Vab 
being negative in this case. The image dipoles A' and B' 
are constructed by appropriate reflection of the dipoles 
A and B. The resulting signs of the interactions Vab 1 
and Vba' between dipoles and image dipoles are nega- 



tive/positive if the respective dipole moments are paral- 
lcl/antiparallcl. Figure 5 shows two electric dipoles near 
a purely electrically polarizable sphere in the triangular 
configuration together with their images in the sphere. 
When the intcr-dipolc angular separation is very small, 
the curvature of the spherical surface can be disregarded 
and the sphere can be approximately replaced by a half 
space as in Fig. 5(a). It is seen that [/W, which is a 
product of one indirect and one direct interaction, is pos- 
itive. Since the negative is a product of two indirect 
interactions, and the direct interaction is stronger than 
indirect one for small interatomic separations the sum 
+ is positive and hence the total potential Uab 
is weakened compared to that in free space. This con- 
firms the numerical results for short distances presented 
in Fig. 3(a). The case of two dipoles located at the oppo- 
site ends of a sphere diameter is sketched in Fig. 5(b). It 
can be seen that f/W is negative. As a consequence, Uab 
is enhanced in agreement with the curves in Fig. 3(a). 

The case of two electric dipoles near a purely mag- 
netically polarizable sphere can be treated by consider- 
ing two magnetic dipoles near a purely electrically po- 
larizable sphere, since the two situations are equivalent 
due to the duality between electric and magnetic fields. 
Again we consider the triangular configuration. In the 
limit of small separation angles, the surface can be re- 
garded as flat. The sketch of the dipoles and their im- 
ages in Fig. 6(a) indicates a negativeness of leading 
to an enhancement of the overall interaction potential [cf. 
Fig. 3(b) inset, solid curve]. When the separation angle 
is large (atoms located on opposite sides of the sphere) , 
it can be inferred from Fig. 6(b) that is positive, 
resulting in a reduction of the interaction potential, in 
agreement with the numerical results shown in Fig. 3(b). 
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FIG. 7: Two dipoles near a sphere are shown in a linear con- 
figuration. Figure 7(a) is showing two electric dipoles while 
in Fig. 7(b) two magnetic dipoles are shown. 



We turn now to the linear configuration. For dipoles 
situated near a purely electrically polarizable sphere 
[Fig. 7(a)], is negative resulting in an enhancement 
of the total interaction potential for all distance regimes 
as visible in Fig. 4(a). For a purely magnetically po- 
larizable sphere, we again invoke the duality principle to 
replace it by a purely electrically polarizable one, and the 
electric dipoles by magnetic ones as shown in Fig. 7(b). 
It can be inferred from the sketch that is positive for 
all distances. In order to be conclusive about the body- 
induced effects, one hence has to compare the magnitudes 
of the competing and U^- 2 \ For small atom-atom 
separations, the direct interaction dominates, so is 
stronger than and the potential is reduced as shown 
in Fig. 4(b), inset. As the interatomic separation in- 
creases, the indirect interaction gains in relevance and 
hence may become dominant leading to an enhance- 
ment of the total vdW potential, in agreement with the 
curves presented in Fig. 4(b). 



V. SUMMARY 

We have studied the mutual vdW interaction between 
two atoms near a dispersing and absorbing magneto- 
electric sphere and presented both analytical and nu- 
merical results. When the radius of the sphere becomes 
sufficiently large, then the interaction potential tends to 
the one found for two atoms near a magneto-electric half 
space. In the opposite case of a very small sphere, the 
sphere can be regarded as being a third atom with respec- 
tive electric and magnetic polarizabilitics. In particular 
for electrically polarizable atoms, the three-atom inter- 
action potential is recovered. 

The numerical calculations performed for medium- 
sized spheres show that — compared to the case of the 
atoms being in free space — the interatomic vdW interac- 



tion can be enhanced as well as reduced, depending on 
the electromagnetical properties of the sphere, the posi- 
tions of the atoms relative to the sphere, and the posi- 
tions of the atoms relative to each other. In general, the 
electric properties of the sphere have a more pronounced 
influence on the potential than the magnetic properties. 
We have shown that the behavior of the interatomic po- 
tential can be qualitytively understood in the basis of an 
image-charge model. 

The results also indicate an essential difference be- 
tween finite- and infinite-sized systems, which particu- 
larly becomes apparent in the linear configuration with 
a purely electrically polarizable sphere. Depending on 
the distance between the surface and the neighbouring 
atom, the (normalized) potential here approach different 
values, while in the case of a half-space they converge to 
a single value. 
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APPENDIX A: DERIVATION OF EQS. (24)-(26) 

To perform the summations over m and p in Eq. (17), 
we begin with the case p=—l and evaluate the first term 
in the square brackets. Using Eqs. (18) and (22), we may 
write 

V C nm M nmi - 1 (r A ,k)M nmi - 1 (r B ,k) = ?" 

^— sin Oa sin B 

m=0 

n 

x °nmm 2 cos(TO7r)P™(cos^)P„ m (cos0 s ), (Al) 

m=0 

where 

C nm = fr*~ m i' (2-(W). (A2) 
(n + my. 

Differentiating the addition theorem for spherical har- 
monics 

n 

J2 C nm P™ (cos 9) P™ (cos 0') cos(m^) = P„(V>), (A3) 
where 

ip = cos cos 0' + sin sin 0' cos 0, (A4) 
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twice with respect to </>, we obtain 

n 

C nm m 2 cos(to0)P™(cos6»)P™(cos0') 



m=0 



d 2 P n {^) (d^\ 2 dP n (V>) d 2 V> 

d</> 2 \d(t>) dtp dcf) 2 ' 1 ' 

Using Eq. (A5) together with Eq. (A4) in Eq. (Al), we 
find that 



m=0 



CnmM„ TO) _i(r A ,A;)M„ m) _i(rB,fc) 

_ _ n (o) dP„( 7 ) 
~ Wn d 7 



e eA e eB (A6) 



(recall that 7 = cos , 6 = 0,4 + #_b ) • The summation over 
m in the other terms in Eq. (17) can be performed in a 
similar way to obtain 



C nm M nm!+1 (r A ,k)M (rs, fc) 

= Ql°)F n ( 7 )e ,e 0B , (A7) 



m=0 



C„ TO N„ m _i(r A ,A;)N„ m _i(r B ,fc) 



m=0 



g(f)p;( 7 )e^ A e^, (A8) 



^ ' CnmNnm, + l( r A!^)N n m,+l 



m=0 



-{Kn + 1)] 2 QWP„( 7 ) 



6r d G r 



uj'r A r B 

+ n(n + 1) sin OP; (7) \Q^e rA eg B + Q^e 0A e r 



Ql 4 )p„( 7 )e eA e eB }. (A9) 



Inserting Eqs. (A6)-(A9) in Eq. (17), we arrive at 
Eqs. (24)-(26). 



APPENDIX B: THE POTENTIAL 
CONTRIBUTIONS l/V , EQ.(37), AND u£ } , 
EQ.(38) 



16ir 3 l 5 rArB 



n=l 



/>00 

x / duue _ ' tl/c Q; y i(«u)as(iu)-B^( 
Jo 



. ,0 

m 



(2) 
n 

(3) 



[/ 2 /(0±ff(0^(B)iA(B)], (B2) 



(1) fyigc ^ (2n+l) 



xa A (iu)a B (iu) 



l(m : ''l 5 E n(n + 1 1 /„ 

n— 1 x ' u 

B M (j „ )0 , I)p , (7) _£^M 

2 , 



[l 2 f(0 cos 9 - g(0r A r B sin 2 9] , (B3) 



16tt 3 Z : 



(gn+1) T 



duu 3 e- lu / c 



x a y 4(iM)as(m) 



u"-r A r B 



(B4) 



fc 2 2 00 

L/( 2 )(r A ,r B )= 3 ^" 2 C r2 2 m(m+l)(2m+l) 



x n(n + l)(2n+l)P m (7)P„( 7 ) / duu 2 

Jo 

x a A (iu)a B (iu)B*(iu)B% (iu)Q$Q£\ (B5) 



32^3,2^2 



J2 (2m+l)(2n+l) 



UW(t a ,t b ) 

Or 

/>0O 

xP L {l) P n{l) \ duu 2 a A (iu)a B {iu) 
Jo 

o (2) o 

Q rn Q 



r A r B 1 

^ D m,n=l 



xB»WBZM%%, (B6) 



We substitute Eqs. (32)-(35) together with Eqs. (24)- 
(26) into Eqs. (37) and (38) and obtain the following 

expressions for the nonzero U-p and u[p: 



u£>(t a ,t b ) 



lQir 3 l 5 r A r B 



^n(n+l)(2n+l)P„( 7 ) 



n=l 



/>oo 

/ duu e~ lu/c a A (iu)a B (iu)B% (iu)Q^ 
Jo 

x [/*/(£) cos Q-g(0 UIb], (Bl) 



U^{v A ,v B ) 



hi4 ^ (2n+l) (2m+l) 
n(n+l) m(m+l) 



327r 3 c 



E 



/•OO 

/ duw ( 
Jo 



a^(iw)a B (iu) 



c 2 ^ w P ro ( 7 ) 

u 2 r A r B ^» p 4(7 ) 



J}Mq(\) Pm(l) 
P m ( 7 ) 

" y " F„( 7 ) 
c 2 P^ nW J-n(7) 



u 2 r A r B ^ n P> n (j) 



(B7) 
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APPENDIX C: THE LIMITING CASES OF A 
LARGE AND A SMALL SPHERE 

When in the case of a large sphere the conditions (39) 
and (40) are satisfied, then the leading contributions to 
the sums in Eqs. (24)-(26) come from terms with n>l 
(also see Ref. [31]), for which the spherical Bessel and 
Hankel functions approximate to [35] 



where 



jn{z) 



1 



(2n + l)!! 
and 



+ ■ 



4n + 6 (16n + 24)(2n + 5) 

(CI) 



1 + - 



; + - 



4n-2 (16n-8)(2n-3) 

(C2) 

respectively. Hence, Eqs. (20) and (21) approximate to 



(iu) = 

and 
B»(iu) = 



i(iRu/c) 



2n+l 



16n[(2n+ l)!!] 2 [/x(iu) + 1] 

i(iRu/c) 2n+1 
16n[(2n+ l)!!] 2 [e(m) + 1] 

respectively, where 



(a 2 n 2 +ain + a ) 
(C3) 



(b 2 n 2 + b 1 n+b ) , 
(C4) 



a = {488 + [1 + e(iu)n(iu)] 2 R 4 u 4 /c 4 } [n(iu) - 1] 

- 8 {7fj,(iu)[e(iu)-l]-5e(iu)fj,(iu) 2 +5} R 2 u 2 /c 2 , 

(C5) 

ai = 8 {26 + [1 + e{iu)n(iu)]R 2 u 2 /c 2 } \fj,(iu) - 1] 
= Sr/(iu)\p(iu) — 1], (C6) 



a 2 = 32[/z(iu) - 1], 



(C7) 



and bo, bi, and b 2 can be found from a , Oi, and a 2 , 
respectively, by interchanging fi and e. Equations (27)- 
(30) then approximate to 



<tf> = 



ic\ 2n+2 [(2n-l)!!] 2 
u 



(C8) 



{rATB) n+1 

Q^=Q^ = -Q^/n = -nQ^. (C9) 

In order to illustrate the application of the approxi- 
mation scheme to the Green tensor elements (24)-(26), 

let us consider the element G^J (ta, *b, iu). Inserting 
Eqs. (C4) and (C8) in Eq. (24), we find 



G$(rA,r B ,Mi) = 



-cH 



128iru 2 R 3 {e(iu) + 1] 



^ afcSfc, 

(CIO) 



*k = $> fc t n+1 P„(7) 



71=1 



and t = R 2 1 {jaTb)- Differentiating the identity 



^1 - 2i 7 + i 2 



(Cll) 



(C12) 



[35] with respect to t, we can perform the summation in 
Eq. (CIO) to obtain 



so 



sjl - 2t 7 + t 2 



-t, 



si = 



t 2 j - t 3 



S2 = 



[1 - 2t 7 + i 2 ] 3/2 

3 7 2 - A-yt + 2t 2 - 1 
[1 - 2i 7 + t 2 } 



5/2 



(C13) 
(C14) 
(C15) 



Recalling the conditions (39) and (40), we can further 
simplify the result. Up to second order in the small pa- 
rameters Sa'/R, we have 

* = 1 _ k ^ + W + V*A + *B + ^ (C16) 
implying that 



i-2t 7 + ^e 2 + ( ^ + ^ )2 12 



R 2 



Using Eqs. (C13)-(C17) in Eq. (CIO), we find that within 
this order, 



GW(r A ,r B ,iu) 



1287ra 2 ^[e(m) + 1] 



r si 2 

x I 32{X 2 - 28 2 + )[e{iu) - 1] - %-^±r)(iu)[e(iu) - 1] 



488+ [1 + e(iu)n{iu)} 2 



c 4 



7e(iu) [n(iu) — 1] — 5s 2 (iu)fj,(iu) + 5 



- 1] 



i? 2 u 2 " 


ill 


c 2 


i? 2 J 



(C18) 



Recalling that X, l + , S + <C R, it can be seen that un- 
less \e(iu) — 1| <C 1, the third term in the curly bracket in 
Eq. (C18) can be approximately ignored. Hence, 

G n(, A ,^, w ) = _^__ 



16i? 



[26 + s{iu)n(iu)R 2 u 2 /c 2 ^ 



(C19) 
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In the case of a purely magnetically polarizable sphere 
[e(iu) = 1], the leading contribution to Grr comes from 
the third term in the curly brackets in Eq. (C18): 



G?W(rA,rB,ju) 



jj,(iu) 



16ttZ + 



(C20) 



The other Green tensor elements can be evaluated in a 
quite similar way. Substituting the resulting expressions 
into Eqs. (37) and (38), and summing them in accordance 
with Eq. (36), we eventually arrive at Eqs. (41) and (43). 

In the limiting case of a small sphere where the condi- 
tion (44) holds, the leading contributions to the frequency 
integrals in Eqs. (B1)-(B7) come from the region where 
ii<c/ft, or equivalently \zq\, \zi \ Cl (also see Ref. [31]). 
In this region we may approximate the spherical Besscl 
and Hankcl functions appearing in Eqs. (20) and (21) by 
their ncxt-to-lcading order expansions in z [35], i.e., 



and 



3n{z) 



/#>(*) = - 



(2n+l)!! 
(2n- 1)!! 

7 n+l 



1 



4n + 6 



2-4n 



(C21) 



(C22) 



so that Eqs. (20) and (21), respectively, approximate to 



revealing that in Eq. (36) the U-p terms are small in 

comparison to the uj^ terms and can be neglected, so 
that, in leading order of Ru/c, 



£ Ug\r A ,r B ). 



(C25) 



i,j—r,8,(f> 



Further, it can be seen that in the sums in Eqs. (24)-(26) 
the terms with n = 1 are the leading ones, for which 



(1) _ (l + o)(l + 6) b 



(2) _ (l + o)/(6) „ 

«1 ~ a 2 b 2 



(3) _ f(a)(l + b) b 



QX 



(4) _ f(a)f(b) b 



a 2 b 2 



(C26) 



(C27) 



(C28) 



(C29) 



B™(iu)=i 



2n+ 1 



[(2n + iy.\] 2 [fi(iu)n + n + 1] 

. 2n+l 



2n+3 



(C23) 



and 

B%(iu)=i 



2n+l 



[(2n+ l)U] 2 [e(iu)n + n + 1] 



x{[eM-l](n + l) (^f +1 + 0(^) 2 " +3 }, 



(C24) 



F 1 ( 7 ) = P 1 ( 7 )=7, 



>M/^a _ 2 A*(* u ) - I (Ru 



B?{iu) 



3 /x(iu) + 2 \ c 



2 e(m) - 1 f Ru\ 



3 e(m) + 2 \ c / 



(C30) 



(C31) 



(C32) 



Substituting Eqs. (B1)-(B4) together with Eqs. (C26)- 
(C32) in Eq. (C25), we arrive at Eq. (45). 
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